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•  Using  advanced  property  models,  errors  of  common  ideal  assumptions  are  quantified. 

•  Energy  analyses  study  pressurized  electrolysis  versus  product  gas  compression. 
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We  report  on  a  thermodynamic  analysis  for  water  electrolysis  from  normal  conditions  (P  =  0.1  MPa, 
T  =  298  K)  up  to  heretofore  unaddressed  temperatures  of  1000  K  and  pressures  of  100  MPa.  Thermo¬ 
neutral  and  reversible  potentials  are  determined  using  equations-of-state  published  by  the  International 
Association  for  the  Properties  of  Water  and  Steam  and  the  National  Institute  of  Standards  and  Tech¬ 
nology.  The  need  for  using  accurate  property  models  at  these  elevated  temperatures  and  pressures  is 
exemplified  by  contrasting  results  with  those  obtained  via  ideal  assumptions.  The  utility  of  our  results  is 
demonstrated  by  their  application  in  an  analysis  comparing  pressurized  electrolysis  versus  mechanical 
gas  compression.  Within  the  limits  of  our  analysis,  pressurized  electrolysis  demonstrates  lower  energy 
requirements  albeit  with  electrical  work  composing  a  greater  proportion  of  the  total  energy  input. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Hydrogen  is  expected  to  play  a  critical  role  in  future  energy 
systems,  by  complementing  or  replacing  fossil  fuels  in  numerous 
applications  [1,2].  Beyond  its  application  in  electrochemical  and 
combustion  systems  as  a  zero-emissions  fuel,  hydrogen  also  serves 
as  a  feedstock  in  many  industrial  chemical  processes  (e.g.,  the 
Fischer-Tropsch  processes  for  liquid  hydrocarbon  synthesis.) 
Currently,  the  largest  portion  of  global  hydrogen  production  in¬ 
volves  the  reformation  of  fossil  fuels  (mostly  natural  gas).  Among 
the  alternative,  clean,  hydrogen  pathways  is  water  electrolysis 
powered  by  renewable  or  other  carbon-free  sources. 

We  report  on  a  thermodynamic  analysis  for  water  electrolysis 
from  normal  conditions  ( p  =  0.1  MPa,  T  =  298  K)  up  to  1000  K  and 
100  MPa.  The  reversible  and  thermoneutral  cell  potentials  as 
functions  of  temperature  and  pressure  shall  be  described,  and  their 
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utility  demonstrated  via  energy  analyses  of  archetype  systems 
featuring  competing  pressurization  strategies. 

Operation  outside  normal  conditions  is  foreseen  to  be  advan¬ 
tageous  in  terms  of  total  energy  consumption  and  energy  utiliza¬ 
tion.  Elevated  temperature  decreases  the  reversible  work  required 
to  perform  electrolysis,  making  cogeneration  in  thermal  power 
plants  viable  [3j.  Elevated  pressure  eliminates  the  need  to 
compress  product  gases,  at  the  expense  of  requiring  a  pressurized 
feedwater  input.  Hence,  for  an  electrolysis  plant,  the  lower  energy 
required  for  pressurizing  the  denser  water  may  improve  overall 
efficiency.  Additionally  are  possible  reductions  in  electrochemical 
overpotentials;  improved  reaction  kinetics,  and  enhanced  mass 
transport  performance  especially  in  the  supercritical  phase.  These 
overpotentials  are  however  beyond  the  scope  of  this  work. 

Contemporary  systems  for  hydrogen  production  via  electrolysis 
can  be  classified,  like  their  fuel  cell  counterparts,  according  to  the 
ionic  conductor  type  and  the  temperature  of  operation.  Proton- 
exchange  membrane  electrolysis  cells  (PEMECs)  operate  between 
70  and  100  °C,  alkali  cells  80—90  °C,  and  solid-oxide  electrolysis 
cells  (SOECs)  at  higher  temperatures  (600—1200  °C)  [4—6];  all  of 
which  have  operating  pressures  near  atmospheric.  Experimental 
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efforts  have  explored  the  performance  outside  these  conventional 
envelopes.  For  example:  steam  temperatures  in  proton-exchange 
membrane  (PEM)  systems  using  alternative  membranes  [7];  and 
pressurized  PEMEC  and  SOEC  systems,  having  achieved  pressures 
of  7  MPa  and  1  MPa  respectively  [5,6],  Alkaline  systems, 
approaching  10]  and  in  one  case  exceeding  supercritical  water 
conditions  11]  have  been  reported. 

Theoretical  treatments  that  are  complementary  to  experimental 
efforts  are  available.  In  an  attempt  to  remedy  perceived  inaccura¬ 
cies  in  the  literature  of  the  time,  Leroy  et  al.  calculated  theoretical 
potentials  for  alkaline  electrolysis  in  the  ranges  of  25—250  °C  and 
1—100  atm  [12].  The  expressions  for  enthalpies  considered  real 
gases,  but  ideal  mixtures;  whereas  expressions  for  Gibbs  free  en¬ 
ergy  (and  thus  reversible  potential)  treated  all  species  as  ideal.  The 
implications  of  different  systems  level  operating  strategies  were 
not  addressed.  Onda  et  al.  calculated  theoretical  potentials  for 
(non-specific)  electrolysis  in  the  ranges  of  25—250  °C  and 
1—700  atm  [13];  applying  these  to  high-pressure  electrolysis  sys¬ 
tems  analysis.  Their  work  followed  the  expressions  developed  by 
Leroy,  though  with  a  different  source  for  reference  pressure  prop¬ 
erties  as  a  function  of  temperature. 

In  the  present  study  we  address  current  challenges  and 
knowledge  gaps  of  the  literature  by:  expanding  the  parameter 
ranges  considered,  incorporating  advanced  property  models, 
demonstrating  the  error  when  inappropriate  models  are  imple¬ 
mented,  and  analysing  electrolysis  configurations  operating  at  the 
extremes  of  the  ranges. 

2.  Method 

The  treatment  assumes  a  reversibly  operated  electrolysis  cell 
(i.e.  without  overpotentials);  where  the  inlet  stream  is  pure  water 
of  a  single  phase  at  prescribed  pressure  and  temperature.  The  two 
independent  outlet  streams  are  pure  hydrogen  and  oxygen  at  the 
same  pressure  and  temperature.  This  is  shown  schematically  in 
Fig.  1.  The  full  and  half-cell  reactions  for  an  (acid  electrolyte)  elec¬ 
trolysis  cell  are: 

H20— ►0.5O2  +  2e~  +  2H+  anode  (1) 

2H+  +  2e~— >H2  cathode  (2) 

H20— >H2  +  0.5O2  cell  (3) 

The  net  energy  input  (work  &  heat)  required  is  given  by  the  1st 
law  of  thermodynamics.  Neglecting  potential,  kinetic  and  other 
terms,  it  may  be  expressed  for  a  control  volume  (CV)  encompassing 
the  cell  operating  at  steady-state: 

d£cv  ... 

-g—  =  0  =  Q-W+  nH2ofr(H2o|p,r)  -  ~  |io2fyo2|  p,T) 

(4) 
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Fig.  1.  System  schematic. 


where  h^p  T  denotes  molar  enthalpy  at  temperature  T  and  pressure 
p  for  species  i.  Rearranging  and  scaling  terms  per  coulomb  of  charge 
transferred  yields: 


_  W  -  Q  _  h(H20|p,r)  ~  ft(H2|P,T)  ~  0-5fr(o2|p,r) 

tlh  iiH2nF  nF 

where  £th  is  the  thermoneutral  cell  potential. 

The  reversible  work  input  required  is  obtained  by  substituting 
the  2nd  law  of  thermodynamics  for  the  Q  term  of  Eq.  (4).  What 
follows  is  an  expression  involving  the  Gibbs  free  energies  of  the 
species: 


w  =  hH2o (h(H20|p,T)  “  TS(H2OIp,T))  -  4h2  (h(H2|p,T)  -  7s(H2|p,T)) 

~no2  (fyo2| p,T)  -  Ts( o2|p,d)  (6) 

Rewriting  the  equation  per  coulomb  of  charge  transferred  yields 
the  reversible  cell  potential. 


C  W  g(H20|p,T)  -£(H2|p,T)  -  °-5^(02jp,r) 

£rev'cel>  =  rVtF  = - nF -  (7) 

The  quality  of  a  (thermo  analysis)  model  is  contingent  on  the 
accuracy  of  supporting  data.  This  work  makes  use  of  the  most 
recent  equations-of-state  (EOS)  for  thermodynamic  properties 
recognized  by:  the  International  Association  for  the  Properties  of 
Water  and  Steam  (IAPWS),  the  National  Institute  of  Standards  and 
Technology  (NIST),  and  the  International  Council  for  Science  : 
Committee  on  Data  for  Science  and  Technology  (CODATA).  More 
specifically,  the  following  sources  are  used: 

•  For  water,  this  analysis  uses  the  IAPWS  1995  release  on  the 
properties  of  ordinary  water,  as  implemented  in  the  WATER95 
MATLAB  library  [14,15],  IAPWS-95  describes  water  over  the 
ranges  of  the  melting  pressure  curve  to  1273  K  at  1000  MPa, 
though  this  may  be  extended  to  100  GPa  and  5000  K. 

•  For  (normal)  hydrogen,  this  analysis  uses  an  expression  devel¬ 
oped  by  Leachman  et  al.,  valid  up  to  1000  K  and  2000  MPa  [16], 

•  For  oxygen,  this  analysis  uses  an  expression  developed  by 
Schmidt  and  Wagner  [17],  Originally  conceived  for  an  upper 
limit  of  300  I<  and  818  bar,  it  may  be  extrapolated  with  accuracy 
to  1000  I<  and  100  MPa  [18], 

•  For  translation  of  the  aforementioned  data  to  a  consistent 
reference  frame,  this  analysis  uses  standard-state  enthalpy  of 
formation  and  entropy  from  CODATA  Key  Values  for  Thermo¬ 
dynamics  [19], 

The  aforementioned  equations-of-state  describe  succinctly 
the  thermodynamic  properties  of  their  respective  fluids  via  two- 
part  expressions  for  non-dimensional  Helmholtz  free  energy  as  a 
function  of  reduced  density  and  reduced  inverse  temperature 
(Eq.  (8)). 


<P(5,t)  =  cp°(5,T)  +  </>r(5,T)  where  :  c P  =  ^  6  =  y  T  =  y 

(8) 

Such  description  is  advantageous  as  it  does  not  require  ancillary 
equations  (e.g.:  is  not  susceptible  to  errors  from  multiple  equation 
incongruities).  Appropriate  partial  derivatives  of  the  ideal  and  re¬ 
sidual  Helmholtz  free  energy  components  yield  the  necessary  pa¬ 
rameters  (Eqs.  (9)— (11)). 
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p  =  pRT 


1  +<5 


35 


(9) 


The  gas  and  water  compressors  are  taken  as  adiabatic  and 
isentropic.  Eq.  (14)  expresses  the  work  flux  for  these  devices. 


h_ 

RT 


(10) 


W 


compressor 
1  — >2 


=  ni(h(i\p1,T),-h(i\p2j2)j,  where :s(i|PliTl)=s(i|p2ir2) 

(14) 


To  emphasize  the  importance  of  considering  real  thermody¬ 
namic  properties,  reversible  and  thermoneutral  cell  potentials  are 
also  calculated  using  ideal  (but  non-perfect)  assumptions.  For  this 
scenario,  if  the  water  is  in  liquid  phase,  real  IAPWS  data  are  used. 
Otherwise  water,  hydrogen  and  oxygen  enthalpies  and  Gibbs  free 
energies  are  calculated  per  Eqs.  (12)  and  (13): 

h(pJ)  =  h(pre{,T )  (12) 

g(p,T)  =  h(pref,r)  -  Ts(p,T) 


where  enthalpy  and  entropy  as  functions  of  temperature  at  refer¬ 
ence  pressure  are  calculated  from  NIST/IAPWS  data;  whilst  CODATA 
data  are  again  used  for  setting  the  species  independent  reference 
frame. 

The  second  part  of  this  work  is  a  limited  investigation  of  two 
illustrative  strategies  for  obtaining  product  gas  in  a  compressed 
state.  The  first  strategy  (A)  operates  an  electrolysis  cell  in  a  pres¬ 
surized  state  by  running  feedwater  through  a  pump.  The  second  (B) 
operates  an  electrolysis  cell  at  atmospheric  pressure  with  subse¬ 
quent  mechanical  compression  of  product  gases  (including  the 
oxygen).  Provisions  are  allowed  to  accommodate  elevated  oper¬ 
ating  temperature.  The  two  configurations  are  presented  in  Fig.  2. 


The  heat  exchange  elements  regulate  temperatures  prior  to  the 
electrolyser  and  cool  product  gases.  To  reduce  complexity;  heating 
is  taken  as  an  energy  requirement  into  the  system,  whereas  cooling 
occurs  freely  to  the  surroundings  (i.e.:  no  economizers  or  similar 
thermal  energy  recovery  systems).  Eq.  (15)  expresses  the  heat  flux 
for  these  devices. 

Qw2  =  rii(*(i|p,,r2)  -  tyi|p,,T,))  (15) 

Work  and  total  energy  potentials  are  calculated  for  each  strat¬ 
egy.  These  are  the  energy  fluxes  scaled  to  per  coulomb  charge 
transferred.  The  work  comprises  contributions  from  the  compres¬ 
sors  and  the  electrolysis  reaction.  The  total  energy  includes  this 
work  and  contributions  from  external  heat  and  the  endothermic 
component  of  electrolysis.  These  calculated  following  Eqs. 
(16)— (19). 

F  S( h20|A3)  -  £(h2|A4)  ~  °-5g(o2|A4) 

c  work  np 

strategy  A 

,  fythOjA,)  ~fr(H20|A2) 

nF  1 

f  V20|A3)  -  h(H2|A4)  -  °-5h(02!A4) 

L  total  np 

strategy  A 

h(H2Q\A,)  -  ft(H2Q|A2)  V2Q|A3)  -  ft(H2Q|A2) 

nF  nF 

(17) 


(A) 


Qh 


eater  A 


Q 


Q  Q 


(B) 


Fig.  2.  System  schematics.  Above:  pressurized  electrolysis  system.  Below:  ambient  pressure  electrolysis  with  adiabatic  mechanical  product  compression. 
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F  =g(H20|B2)  ~g(H2|B3)  ~  0-5g(O2[B3) 

L  work  np 

stragegy  B 

h(H2|B4)  ~  V2|B5)  ft(Q2|B4)  ~  h(02|B5) 

nF  2  nF 

c  _  h(H20|B2)  -  h(H2|B3)  -  0  5h(02|B3) 

c  total  np 

strategy  B 

/l(H2|B4)  ~  ft(H2|B5)  h(Q2|B4)  -  ^(Q2|B5) 

nF  2  nF 

^(H2Q|B2)  -  ftcHaOlB,) 

nF 


(18) 


(19) 


Fig.  4.  Thermoneutral  and  reversible  cell  potentials  plotted  vs.  pressure  for  various 
temperatures. 


3.  Results  and  discussion 

The  thermoneutral  and  reversible  cell  potentials  are  plotted  as  a 
function  of  temperature  for  several  pressures  in  Fig.  3.  Analogously, 
thermoneutral  and  reversible  cell  potentials  are  plotted  as  a  func¬ 
tion  of  pressure  for  several  temperatures  in  Fig.  4.  Note  the 
conformance  to  IUPAC  conventions  (where  the  potential  of  an 
electrolytic  cell  is  negative). 

As  temperature  increases,  the  reversible  potential  decreases  in 
magnitude.  The  required  work  is  being  offset  by  the  thermal  energy 
contribution  to  the  reaction;  which  is  the  difference  of  the  revers¬ 
ible  and  thermoneutral  potentials,  and  generally  increases  with 
temperature.  The  trend  of  the  total  energy  requirement  (thermo¬ 
neutral  potential)  with  temperature  is  pressure  dependent;  with 
lower  pressures  yielding  increases  in  magnitude  as  temperatures 
increases,  and  the  inverse  for  higher  pressures.  The  discontinuities 
in  thermoneutral  potential,  attributable  to  feedwater  being  deliv¬ 
ered  as  either  liquid  or  gaseous,  are  present  because  this  analysis 
does  not  include  the  process  steps  of  taking  water  (at  arbitrary 
conditions)  to  specified  pressure  and  temperature. 

As  system  pressure  increases,  the  reversible  potential  increases 
in  magnitude;  most  pronounced  for  the  first  few  pressure  in¬ 
crements.  The  change  in  thermoneutral  potential  as  pressure  in¬ 
creases  is  less  intuitive.  At  lower  temperatures  (e.g.:  298  K),  the 
magnitude  first  decreases  then  increases;  whereas  at  higher  tem¬ 
peratures,  the  thermoneutral  potential  increases  in  magnitude 
monotonically. 

Comparing  reversible  and  thermoneutral  potentials;  we  observe 
a  substantial  increase  in  the  work  input  required  for  electrolysis  (to 
jump  say  99  bar  from  1  bar),  whilst  the  net  energy  requirement 


Fig.  5.  Reversible  work  as  a  percentage  of  total  energy  into  cell  vs.  pressure. 

changes  only  marginally.  This  suggests  a  shift  occurs  in  the  distri¬ 
bution  of  energy  input  requirements,  from  thermal  (heat)  to  elec¬ 
trical  (work).  This  phenomenon  is  illustrated  in  Fig.  5,  which  plots 
the  ratio  of  reversible  and  thermoneutral  potentials  (equivalent  to 
the  percent  work  in  the  net  energy  for  electrolysis).  Such  behaviour 
is  not  however  consistent  over  the  entire  analysis  range.  At  higher 
pressures  and  temperatures  (say  increasing  from  700  to  900  bar  at 
800  K);  the  increase  in  work  needed  is  smaller  and  similar  in 
magnitude  to  the  change  in  total  energy  necessary.  Here,  the 
electrical  versus  thermal  balance  is  shifting  back  towards  thermal. 
Had  ideal  gas  assumptions  been  made,  as  is  prevalent  in  the 


Fig.  3.  Thermoneutral  and  reversible  cell  potentials  plotted  vs.  temperature  for  various 
pressures. 


Temperature  /  K 


Fig.  6.  Contour  plot  of  percent  error  in  thermoneutral  cell  potential  vs.  pressure  and 
temperature. 
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Temperature  /  K 

Fig.  7.  Contour  plot  of  percent  error  in  reversible  cell  potential  vs.  pressure  and 
temperature. 

literature,  enthalpy  being  a  function  of  temperature  alone  would 
have  failed  to  reveal  these  behaviours. 

Finally,  to  emphasise  the  need  to  consider  real  species  proper¬ 
ties  when  performing  thermodynamic  calculations  at  elevated 
pressures  and  temperatures;  contour  plots  of  the  percent  error  in 
thermoneutral  and  reversible  cell  potentials  with  ideal  assump¬ 
tions  for  gas  species  are  presented  in  Figs.  6  and  7.  The  error  in  the 
thermoneutral  and  reversible  cell  potentials  with  ideal  assump¬ 
tions  for  gaseous  species  are  greater  than  10%  and  22%  over  the 
investigated  range  of  pressure  and  temperature.  Under  conditions 


where  water  is  in  a  liquid  state,  the  error  is  less  as  real  water  data 
are  used.  Were  the  additional  simplification  of  perfect  gas  (i.e.: 
constant  specific  heats)  behaviour  been  made,  even  greater  error 
would  be  expected  as  enthalpy  and  entropy  would  not  track 
accurately  with  temperature.  Incorrect  determination  of  thermo¬ 
neutral  and  reversible  potentials  can  impact  directly  cell  perfor¬ 
mance  predictions;  at  the  limits  of  our  temperature  and  pressure 
(and  ignoring  irreversibilities),  the  22%  error  in  reversible  potential 
would  yield  an  approximate  7.2  kW  h  overestimate  in  the  elec¬ 
tricity  required  to  produce  1  kg  of  hydrogen. 

Energy  requirements  for  strategies  A  and  B  for  product  gas 
synthesis  and  compression  are  presented.  These  are  reported  as 
total  energy  and  work  potentials,  which  are  analogous  to  reversible 
and  thermoneutral  potentials  of  the  electrolysis  cell;  but  with  the 
inclusion  of  fluid  compression  and  heating  contributions.  These 
potentials  are  plotted  in  Fig.  8. 

Both  the  total  energy  and  work  requirements  are  consistently 
lower  (at  a  given  temperature  and  storage  pressure)  for  pressurized 
electrolysis  (A)  versus  ambient  pressure  electrolysis  (B)  followed  by 
mechanical  gas  compression.  The  total  energy  premium  for  oper¬ 
ating  strategy  A  at  progressively  higher  pressures  is  miniscule 
compared  to  that  of  strategy  B.  In  fact,  studying  the  sub-critical 
temperature  curves  of  298  or  400  I<  from  (a)  in  Fig.  8,  a  minimum 
is  observed  as  pressure  increases.  The  total  energy  for  B  as  a  function 
of  pressure  is  defined  exclusively  by  increasing  compressor  work,  as 
both  heating  and  electrolysis  occurs  at  reference  pressure. 

Examining  the  percent  ratio  of  work  and  total  potentials,  illus¬ 
trated  in  Fig.  9,  reveals  strategy  A  consistently  requires  a  greater 
fraction  of  its  energy  input  as  work  versus  strategy  B.  Though  the 
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Fig.  8.  Work  and  total  energy  potentials  for  competitive  production  strategies  A  and  B.  Negative  potentials  increasing  in  magnitude  represent  greater  energies  delivered  into  a 
system,  (a)  Total  as  a  function  of  pressure  for  select  temperatures,  (b)  Total  as  a  function  of  temperature  for  select  pressures,  (c)  Work  as  a  function  of  pressure  for  select  tem¬ 
peratures.  (d)  Work  as  a  function  of  temperature  for  select  pressures. 
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Fig.  9.  Percent  ratio  of  work  and  total  energy  potentials  for  production  strategies  A  and 
B  as  a  function  of  pressure. 

integration  with  thermal  facilities  is  an  associated  benefit  of  high 
temperature  electrolysis  (as  suggested  by  the  decreasing  ratios 
with  temperature);  its  viability  should  also  consider  electrolysis 
pressure.  The  difference  between  strategy  A  and  B,  and  therefore 
the  need  to  consider  operating  pressure,  is  more  pronounced  at 
lower  temperatures;  where  an  electrolysis  system  might  otherwise 
be  considered  fit  for  a  bottoming  cycle  in  a  thermal  power  plant. 
This  conclusion  would  only  be  strengthened  for  real  electrolysis 
systems,  because  electrical  energy  expended  into  overpotentials 
and  converted  to  heat  would  further  reduce  the  window  for  inte¬ 
grating  external  heat  sources. 

4.  Conclusion 

The  reversible  and  thermoneutral  cell  potentials  for  water 
electrolysis  were  determined  for  temperatures  and  pressures  up  to 
1000  K  and  100  MPa  respectively.  The  potentials  were  obtained 
with  sophisticated  real  thermodynamic  data  based  on  expressions 
for  species  Helmholtz  free  energy.  For  the  conditions  investigated, 
the  error  of  contemporary  literature  which  makes  use  of  ideal 
thermodynamics  relationships  was  shown  to  be  as  large  as  22%. 

Applying  these  results  in  an  investigation  of  archetype  elec¬ 
trolysis  plants,  featuring  either  pressurized  electrolysis  or  post¬ 
electrolysis  gas  compression,  the  energy  requirements  were 
determined  over  a  heretofore  unaddressed  operating  range  of 
pressure  and  temperature.  Within  the  scope  of  this  study,  pres¬ 
surized  electrolysis  required  consistently  smaller  work  and  total 
(work  and  heat)  energy  inputs.  Further,  the  percentage  of  work 
composing  total  energy  for  pressurized  electrolysis  was  also 
greater;  suggesting  the  prospects  of  integrating  an  electrolysis 
system  with  external  thermal  sources  needs  consider  both  oper¬ 
ating  temperature  and  product  pressurization  strategy  (be  it 


electrochemical  or  mechanical).  Ultimately,  a  design  would  be 
driven  by  the  quantity  and  quality  of  available  energy  sources. 


References 

[1]  L.  Barreto,  A.  Makihira,  K.  Riahi,  Int.  J.  Hydrogen  Energy  28  (2003)  267. 

[2]  C.J.  Winter,  J.  Nitsch,  Hydrogen  as  an  Energy  Carrier:  Technologies,  Systems, 
Economy,  Springer-Verlag  New  York  Inc,  New  York,  NY,  1988. 

[3]  W.  Doenitz,  R.  Schmidberger,  E.  Steinheil,  R.  Streicher,  Int.  J.  Hydrogen  Energy 
5  (1980)  55. 

[4]  S.C.  Singhal,  K.  Kendall,  High  Temperature  Solid  Oxide  Fuel  Cells:  Funda¬ 
mentals,  Design,  and  Applications,  Elsevier  Science  Ltd,  2003. 

[5]  A.  Hauch,  S.D.  Ebbesen,  S.H.  Jensen,  M.  Mogensen,  J.  Mater.  Chem.  18  (2008) 
2331. 

[6]  K.  Zeng,  D.  Zhang,  Prog.  Energy  Combust.  Sci.  36  (2010)  307. 

[7]  D.  Aili,  M.K.  Hansen,  C.  Pan,  Q.  Li,  E.  Christensen,  J.O.  Jensen,  N.J.  Bjerrum,  Int.  J. 
Hydrogen  Energy  36  (201 1 )  6985. 

[10]  J.C.  Ganley,  Int.  J.  Hydrogen  Energy  34  (2009)  3604. 

[11]  H.  Boll,  E.U.  Franck,  H.  Weingartner,  J.  Chem.  Thermodyn.  35  (2003)  625. 

[12]  R.L.  Leroy,  C.T.  Bowen,  D.J.  Eroy,  J.  Electrochem.  Soc.  127  (1980)  1954. 

[13]  K.  Onda,  T.  Kyakuno,  K.  Hattori,  K.  Ito,  J.  Power  Sources  132  (2004)  64. 

[14]  Revised  Release  on  the  IAPWS  Formulation  1995  for  the  Thermodynamic 
Properties  of  Ordinary  Water  Substance  for  General  and  Scientific  Use,  2009. 
Doorwerth,  The  Netherlands. 

[15]  P.  Junglas,  Implementing  the  IAPWS-95  Standard  in  MATLAB®,  2008. 

[16]  J.W.  Leachman,  R.T.  Jacobsen,  S.G.  Penoncello,  E.W.  Lemmon,  J.  Phys.  Chem. 
Ref.  Data  38  (2009)  721. 

[17]  R.  Schmidt,  W.  Wagner,  Fluid  Phase  Equilib.  19  (1985)  175. 

[18]  E.W.  Lemmon,  R.T.  Jacobsen,  S.G.  Penoncello,  D.G.  Friend,  J.  Phys.  Chem. 
Ref.  Data  29  (2000)  331. 

[19]  J.D.  Cox,  D.D.  Wagman,  V.A.  Medvedev,  CODATA  Key  Values  for  Thermody¬ 
namics,  1989.  New  York,  NY. 

Nomenclature 


T:  temperature,  K 
p:  pressure,  Pa 

Tref:  reference  temperature,  298.1  K 
prefi  reference  pressure,  1  x  105  Pa 
p:  density,  mol  m  3 
h;  mole  flux,  mol  s-1 
h:  molar  enthalpy,  J  mol1 
s;  molar  entropy,  J  K  1  mol-1 
g:  molar  Gibbs  free  energy,  J  mol-1 
a:  molar  Helmholtz  free  energy,  J  mol-1 
x//:  molar  exergy,  J  mol  1 
Q:  heat  flux  (inwards  +  ve),  J  s_1 
VV:  work  flux  (outwards  +  ve),  J  s_1 
Erev:  reversible  cell  potential,  V 
Eth:  thermoneutral  cell  potential,  V 
n;  number  of  electrons,  2  e~ 

F:  Faraday  constant,  96485  C  mol  1 

R:  ideal  gas  constant,  (function  of  EOS,  8.314  for  ideal  calculations)  J  K_1  mol-1 

r;  inverse  reduced  temperature,  dimensionless 

cp:  non-dimensional  Helmholtz  free  energy,  dimensionless 

6:  reduced  density,  dimensionless 

Subscript c:  supercritical-point  property,  — 

Superscript0:  ideal  component,  — 

Superscript:  residual  component.  — 


